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CN 1 ABSTRACT 

We present observations and analysis of the broadband afterglow of Swift GRB 071025. 
Using optical and infrared (RIYJHK) photometry, we derive a photometric redshift 
of 4.4 < z < 5.2; at this redshift our simultaneous multicolour observations begin at 
~30 s after the GRB trigger in the host frame, during the initial rising phase of the 
afterglow. We associate the light curve peak at ^580 s in the observer frame with 
the formation of the forward shock, giving an estimate of the initial Lorentz factor 
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■ ro ~ 200. The red spectral energy distribution (even in regions not affected by the 

Lyman-a break) provides secure evidence of a large dust column. However, the inferred 
extinction curve shows a prominent flat component between 2000-3000 A in the rest- 
frame, inconsistent with any locally observed template but well-fit by models of dust 
formed by supernovae. Time-dependent fits to the extinction profile reveal no evidence 
of dust destruction and limit the decrease in the extinction column to AA3000 < 0.54 
mag after t — 50 s in the rest frame. Our observations provide evidence of a transition 
in dust properties at z ~ 5, in agreement with studies of high-z quasars, and suggest 
that SN-formed dust continues to dominate the opacity of typical galaxies at this 
redshift. 

Key words: gamma-rays: bursts — dust, extinction 



1 INTRODUCTION 

Starting with the discovery of the 9th magnitude afterglow 
of gamma-ray burst (GRB) 990123 (|Akerlof et al. 1999), 

* e-mail: dperley@astro.berkeley.edu 



the early-time study of GRB afterglows has presented great 
promise to elucidate both the nature of the gamma-ray 
burst phenomenon itself and of the medium surrounding 
these objects in extremely distant galaxies. Fast-responding 
telescopes, slewing to the burst location in time to catch 
the afterglow at or near the time of peak luminosity, can 
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probe the physics of the explosion in the initial seconds 
as the ultrarelativistic outflow is decelerated by the in- 
terstellar medium. Continued observations can then follow 
the evolution of the reverse and forward shocks for many 
hours as the afterglow fades away, providing constraints 
on the still poorly-understood early-time emission pro- 
cesses . In addition, the extreme luminosities at early times 
(e.g.. iKann et al.1 l2007al ; iBloom et al] 120091 ; iRacusin et ail 
2008) enable even very small telescopes to provide pre- 
cise photometric and occasionally spectroscopic measure- 
ments of the afterglow spectral energy distribution (SED) 
and act as a probe o f interstellar gas and dust out to the 
epoch of reionization llKawai et al.ll2006l; iTotani et al. 20061; 
Gallerani et al.1 120081; iMcQuinn et alj |200S|; iGreiner et al.1 



2009al ; lTanvir et al.ll2009l ; ISalvaterra et al.ll200Sl ). And while 
the usage of early-time SEDs as probes of the interestel- 
lar environment is hindered to some extent by the uncer- 
tain emission processes acting at these times, they never- 
theless can provide constraints on the direct influence of 
a GRB on its surrounding med ium in the form of dust 
destruction and photoionization dWaxman fc Draind 2000: 
Fruchtcr ct al. 200f; Dra ine fc Had 12002 ; iPerna fc Lazzatil 
20021 ; IPerna et , al.ll2003h . 



At the same time, however, the fleeting and time- 
variable nature of GRB afterglows poses several challenges 
for these early-time diagnostics. To maximise sensitivity, 
the smallest telescopes typically do not employ filter sys- 
tems and therefore give minimal frequency-domain infor- 
mation. When filters are employed, ordinary telescopes 
are forced to employ a filter cycle, creating the possi- 
bility of confusion between spectral and temporal evo- 
lution of the event. Nevertheless, progress has advanced 
steadily with the commissioning of several simultaneous- 
colour robotic telescopes. The P eters Automatic In frared 
Imaging Telescope (PAIRITEL; IBloom et ail 120061 ). on- 
line since late 2004, provides simultaneous measurements 
in the J, H, and .Ks-bands every 7.8 s starting within 
1-3 minutes of a typical GRB and is the primary sub- 
ject of this paper. Notable PAIRITEL-followe d bursts in- 
clude GRBs 041219A, 06 1 126, and 080 3 19B (|Blake et al.1 
120051 ; iPerlev et all l2008al ; IBloom et all 12009] ). More re- 
cently, the seven-channel Ga mma-Ray Burst Opt ical/Near- 
Infrared Detector (GROND, fGreiner et aTll2008l) has also 
produced simultaneous SEDs of afterglows at over the 
wavelength range 4000-24000 A, including in several cases 
time-dependen t SED s during the afterglow rise and fall 
Kriihler et"aH |2008|. [2009b) and short-timescale flares 



Kriihler et al 



20081 . 

l2009al ; lGreiner et al.|[2009bh . and RAPTOR- 



T has tracked spectral changes during the fading of 
GRB 080319B in s everal optical bands simultaneously 
<|Wozniak et al.ll2009l ). In all cases, colour evolution appears 
to be absent or modest, consistent with the lack of strong 
colour evolution in the generally less constraining measure- 
ments by filter-cyc ling instruments such as the Swift UVOT 
l|Oates et al.ll2009h . Correlation with the gamma-ray prompt 
emission and with X-ray flares (al so thought to be as- 
sociated with the pr ompt ph ase: [K occvski e t al.1 120071 . 
IChincarini et al.1 120071 ) is rare (lYost et al.ll2007l). but has 



been observed in some cases (IVestrand et al.1 120051 
Page et al. 2007]; iKriihler et al.1 l2009al ; I Racusin et al.l 
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ferent models for the large variety of l ight curve behavior s 
observed at early ti mes: reverse shock (Sari fc Piran|[l999l ). 
energy rei niection (iRees fc Meszarosl 1 19981) . prompt emis- 
sion (e.g^jJKumar & Pan aitescull2000l. 120 08*). outflow decel- 
eration ( Sari fc Piranlll999l ; lMeszarosll2006|). spectral b reaks 



moving through the optical bandpass 1 Sari et alii 19981 ). and 
many others. 

GRB 071025, detected by the Swift mission 
llGehrels et al. I 120041 ). provides among the best probes 
of the early-time behavior of a gamma-ray burst to date. 
While no secure spectroscopic redshift was attained (z ~ 5.2 
was estim ated from a low-qu ality HIRES optical spectrum 
at Keck; iFvnbo et al.l 120091 ). the photometric SED pre- 
sented here shows clear evidence of a Lyman-a break in 
the observer-frame i?-band and indicates a photometric 
redshift of 4.4 < z < 5.2 ( i]3.2[l . making this among the 
highest-redshift bursts to date and one of only a few 
observed in simultaneous colours during prompt emission. 
Our infrared and optical observations start at ~30 s after 
the burst in the rest frame and follow the rise, peak, and fall 
of an afterglow in simultaneous rest-frame optical colours. 
In this paper, we use this unique data set to test various 
models for the origin of the early emission and conclude it 
is likely due to the deceleration of the burst outflow into a 
uniform-density interstellar medium, allowing estimation of 
the Lorentz factor F ( i]4.1[) . The IYJHK spectral energy 
distribution demonstrates the existence of a significant dust 
column obscuring a star-forming region at z > 4.4 and 
provides evidence that the dust at this epoch had different 
properties from dust that prevails along sightlines in the 
more nearby uni verse, in agree m ent w ith the study of a 
high-^ QSO by iMaiolino et all (|2004h . We suggest that 
this difference is reflective of an absence of evolved AGB 
stars in these earliest epochs, and search for (and place 
stringent limits on) signs of destruction of this dust by 
radiation from the GRB ( H4.2p . Throughout the paper we 

and assume cosmological 



use the convention F oc t 



parameters h — 0.71, Qa ~ 0.7, Qt, 



2 OBSERVATIONS 



0.3. 



Klotz et al.1 12009T ) . These multi-band observations are par- 



ticularly important for distinguishing the predictions of dif- 



2.1 Swift 

At 04:08:54 UT on 2007 October 2 5^, the Burst Alert Tele- 
scope (BAT. lBarthelmv et al.l l2005) on-board Swift detected 
GRB 071025 and performed a rapid slew to the GRB loca- 
tion, beginning observations with the XRT (jBurrows et al.l 
120051 ) at 146 s after the trigger. The BAT light curve is broad 
and only slowly variable: the flux rises slowly during the first 
~80 s and peaks several times before beginning a steady de- 
cay at approximately ~130 s. The GRB remains detectable 
above the background until Swift was forced to slew away 
from the position due to an Earth constraint at 422 s af- 
ter the initial trigger. Observations resumed at 3500 s, and 
tracked the afterglow using the XRT with no further large 
gaps in temporal coverage for the next ~3 days, after which 
it became too faint to be detected. Details of ou r high-energy 
reduction pipeline are described in detail by iButler et al.l 



1 This trigger time will be used as the reference time in the re- 
mainder of the paper. 
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|2007l ) for the Swift BAT and bv lButler fc Kocevskil (|2007l ') 
for the Swift XRlfl 

Swift's Ultra , Violet-Optical Telescope (UVOT, 
iRoming et al"1l2005l ) observed the field starting at 155 s, but 
detect ed no significant af terglow signal in any of its seven 
filters (|Pagani et al.ll2007f ). The non-detection is consistent 
with our photometric redshift, as outlined in £13.21 

2.2 PAIRITEL Observations 

The robotic infrared observatory PAIRITEL consists of 
the 1.3-m Peters Telescope at Mt. Hopkins, Arizona, 
formerly used for the Two Micron All Sky Survey 
f2MASS; ISkrutskie et ai1l2006t ), re-outfitted with the south- 
ern 2MASS camera. PAIRITEL, like 2MASS, makes use of 
two dichroics to image in the J, H, and K s filters simulta- 
neously. 

PAIRITEL responded t o the initial Gamma-r ay Burst 
Coordinate Network (GCN, iBarthelmv et al.lll995T ) alert at 
74.3 s and slewed immediately to the source. Observations 
began at 162 s and continued uninterrupted until 3812 s, 
when due to a problem with the observing queue PAIRI- 
TEL temporarily slewed to another location. Observations 
resumed at 9108 s and continued for another two hours. Raw 
data files were processed using standard IR reduction meth- 
ods via PAIRITEL Pipeline III and resampled using SWarp 
l|Bertin et al.ll2002h to create final 1.0-arcsec/pix images for 
final photometry. PAIRITEL 's standard observing cycle is to 
take three 7.8 second exposures in immediate succession at 
each dither position. While the early afterglow is detected in 
even the shortest 7.8 second frames, the S/N was too low for 
reliable photometry, so the shortest exposures reported here 
consist of 23.4-second "triplestacks" , the sum of all three 
images at each dither position. These images were further 
binned at successively later times to further improve the 
S/N. The afterglow position, relative to 2MASS astrometric 
standards, is a = 355.0711583, S = +31.778575 (J2000). 

Photometry was performed in IRAF0 using the phot 
task. Best results were achieved using aperture photometry 
with an aperture radius of 2.25" in J-band, 2.5" in //-band, 
and 2.75" in K 3 -ba,nd. Unfortunately, while conditions dur- 
ing the observations were generally clear, the night was not 
fully photometric, with variations in the transmission of up 
to 0.3 mag during the course of observations and additional 
significant fluctuations in the seeing. Calibration was there- 
fore performed by re-determining the zeropoint for each im- 
age individually by comparison to our secondary field stan- 
dards C CrtTj) . Fortunately, the field of GRB 071025 is rich in 
bright field stars, and a total of eight nearby stars (present 
and well-detected in even short exposures with reference un- 
certainties of <0.05 mag) were used to determine the zero- 
point. The statistical uncertainty on the zero-point (never 
more than 0.05 mag) is essentially negligible relative to other 

2 Swift bursts occurring after these publications, in- 
cluding GRB 071025, have been processed using the 
same methods; these results are available online at 
http: / / astro.bcrkelcy.edu/~nat / swift/ 

a IRAF is distributed by the National Optical Astronomy Ob- 
servatory, which is operated by the Association of Universities 
for Research in Astronomy (AURA) under cooperative agreement 
with the National Science Foundation. 



sources of error in all cases. Systematic sources of error are 
addressed in £13.41 

The large plate scale of PAIRITEL (2.0 "/pix) and the 
variable sub-pixel response function of the NICMOS3 arrays 
creates a significant additional uncertainty in each position 
beyond ordinary photome tric errors, estimated at ~3 per- 
cent bv lBlake et alj (|200&f ). To quantify this uncertainty as 
accurately as possible, we constructed light curves for stan- 
dard stars of different magnitudes in regions of the image 
free of defects by measuring the image-to-image magnitude 
variations of bright (source-dominated) stars. An additional 
uncertainty of approximately ~0.02 mag per position was 
required to incorporate the observed scatter in the photom- 
etry of these objects. Additionally, we examined fainter (sky 
noise-dominated) stars to compare the IRAF-generated un- 
certainty to that observed in the zero-pointed light curve, 
finding the IRAF uncertainties to be too low by about 20% 
in each filter. Therefore the final uncertainties on our pho- 
tometry, reported in Table [3] were determined by increas- 
ing the IRAF uncertainty by 20% and adding the result in 
quadrature with 0.02/^/ 7V pos mag, where N pos is the num- 
ber of unique dither positions per stacked image. 



2.3 REM Observations 

GRB 071025 a lso triggered REM (Rapid-Ey e Mount; 
IZerbi et al.ll200ll) . a robotic (|Covino et al.ll20o4 ) telescope 
located at the ESO Cerro La Silla observatory (Chile). The 
REM telescope has a Ritchey-Chretien configuration with a 
60 cm f /2.2 primary and an overall f /8 focal ratio in a fast 
moving alt-azimuth mount providing two stable Nasmyth 
focal stations. At one of the two foci, the telescope simulta- 
neously feeds, by mea ns of a dichroic, two camer as: REMIR 
dConconi et al.l [20041 ') for the NIR and ROSS (jTosti et al.l 
l2004h for the optical. Both cameras have a field of view of 
10x10 arcmin and imaging capabilities with NIR (1/x, J, H, 
and K) and Johnson- Cousins VRI filters. Observations of 
the GRB 071025 field began at 144 s after the trigger, al- 
though this initial //-band exposure did not detect the after- 
glow. The optical camera was unfortunately not operational 
due to maintenance, so exposures were acquired only in In, 
J, H, and K. 

The raw frames were corrected for dark, bias, and flat 
field following standard procedures. Although the burst was 
at low elevation at the trigger time, seeing conditions were 
good and photometry was performed using a 3.5 pixel (1.2") 
aperture. Conditions were not photometric, and so the zero- 
point was determined for each image individually in JHK 
bands using a subset of 2MASS-based standards. The \\x- 
band (often referred to as z in previous work, though this 
filter has almost no overlap with the traditional SDSS z- 
band), after taking into account the transmission of the 
ROSS/REMIR dichroic, is close to the MKO Y-bamfl and 
so we treated this filter as a Y measurement, using the in- 
terpolated magnitudes in Table[S](see £12. lip and basing the 
calibration on four reference stars well-detected in all im- 
ages. 



ht tp: //www. ukidss.org/tec hnical / instrument / filters.html; see 



alsolfflllcnbran d et ahl (|200: 



4 D. A. Perky et al. 



2.4 RAPTOR Observations 

The RAP TOR (Rapid Telescop es for Optical Response) ex- 
periment l)Vestrand et al] 120021 '). operated by Los Alamos 
National Laboratory, consists of a series of small telescopes 
at the Fenton Hill Observatory in New Mexico. RAPTOR-S 
is a fully autonomous robotic telescope with a 0.4-m aper- 
ture and typical operating focal ratio f/5. It is equipped 
with a 1000x1000 pixel CCD camera employing a back- 
illuminated Marconi CCD47-10 chip with 13/j pixels. 

RAPTOR-S responded automatically to the localiza- 
tion alert and was on target at 04:10:14.95 UT, 81.3 s af- 
ter the trigger time (4.2 s after receiving the GRB posi- 
tion). The rapid response sequence of RAPTOR-S consists 
of nine 5-second images followed by twenty 10-second im- 
ages and finally 170 30-second images for a total of ~2 hours 
of coverage (including 5-second intervals between exposures 
used primarily for readout). In order to improve the S/N 
ratio, photometry was performed on coadded images. Aper- 
ture photometry was p erformed using the SExtractor pack- 
age |Bertin et alj|2002ft . and the magnitude offsets between 
epochs were derived using several dozen field stars. 

Because of the extreme redness of this afterglow, the 
unfiltered RAPTOR-S observations required a special cali- 
bration procedure. Although the effective wavelength of the 
response curve for RAPTOR-S is close to that of the stan- 
dard R band, the spectral energy distribution of this burst 
( i|3.2[) indicates a sharp drop in the flux between I and R 
bands, likely due to the onset of the Lyman-a forest. As a 
result, most photons detected by RAPTOR-S actually fall 
in the spectral region covered by the standard I filter. 

Therefore, we tie the unfiltered data to /-band stan- 
dards from the Lick calibration ( q2.1ip . The offset (mc — 
J)stat between the unfiltered magnitudes and standard I 
was derived using 7 well-measured stars in the vicinity of 
the GRB covering a narrow range of colours 0.5 < (R — I) < 
0.66. Assuming that the SED of the burst emission did not 
change significantly between the time of RAPTOR-S obser- 
vations and the time when it was measured, we derived an 
approximate correction to (mc — /) 8 tar to account for the 
extremely red colour of the G RB. We used a K5V model 
spectrum from iKurucj (|l979T ) as a proxy SED matching 
the mean colour of our comparison stars. By folding both 
SEDs with response curves of RAPTOR-S and the standard 
/-band filter we find (mc — /)grb = (mc — /) 8 tar + 0.74 
mag. The uncertainty of the derived zero point is about 
10%; consistent with this, we measure a relatively small off- 
set of —0.08 magnitudes between the calibrated RAPTOR 
magnitudes and an extrapolation from later, filtered /-band 
observations using our light curve model (see i)3.1l) . Table [4] 
lists the final RAPTOR-S photometry. 



2.5 Super-LOTIS Observations 

Super-LOTIS (Livermore Optical Transient Imaging Sys- 
tem) is a robotic 0.6-m telescope dedicated t o the search for 
optical cou nterparts of gamma-ray bursts (Williams et al. 
2004, 2008). The telescope is housed in a roll-off-roof facil- 
ity at the Steward Observatory Kitt Peak site near Tucson, 
Arizona. Super-LOTIS triggered on GRB 071025 and be- 
gan observations at 04:10:29 UT (95 s after the trigger), 
acquiring a series of _R-band frames, which were reduced 



using standard methods. Unfortunately, because of the op- 
tical faintness of the afterglow and high sky background, the 
quality of the images is poor and even after extensive stack- 
ing the detection is marginal, particularly in the earliest few 
stacks. Photometry was performed using aperture photom- 
etry and our Lick /?-band field calibration as detailed in 

2.6 Lick Infrared Observations 

We acquired an additional series of infrared observations 
using the 3m Shane telescope at Lick Observatory equipped 
with the UCLA GEMINI IR camera (|McLean et all Il993l . 
1 1994) . A total of nine exposures were acquired in J and K' 
bands simultaneously starting at 04:52:23 UT, integrating 
for 1 1 coadds of 20 s each in J and 35 coadds of 6 s each in 
K' . The IR afterglow was still very bright at this time, and is 
clearly detected with signal-to-noise S/N > 50 in individual 
exposures. Reduction was performed via direct subtraction 
of temporally adjacent exposures followed by division by a 
twilight flat. Photometry was performed using IRAF and an 
aperture of 2 pixels (1.4"); images were calibrated relative to 
the PAIRITEL magnitudes of five nearby bright field stars. 

The response of the K' filter is significantly differ- 
ent from K s and the GRB exhibits an apparent colour 
(H — K rs 1.0 that is much redder than any field star used 
for comparisons (ranging between H — K — 0.04 — 0.18). To 
correct to K s for direct comparison to the PAIRITEL data, 
we use an approximate cor rection of K s m K w K' — 0.07 
ijWainscoat fc Cowiel Il992h . with this correction inferred 
from the reddest star in Table 1 of that work (Oph SI, 
H-K = 0.94, K' -K = 0.07±o.oll)- The K to K s colour 
term is assumed to be negligible. This is found to produce 
good agreement between Lick data and coeval PAIRITEL 
points. However, due to uncertain differences between the 
Lick, MKO, and other filter sets and the intrinsic GRB spec- 
trum itself the overall calibration offset could be as much as 
0.05 mag, and as a result the Lick K photometry is not used 
in fitting. 

2.7 MAGNUM Observations 

The MAGNUM (Multicolour Active Galactic NUclei Moni- 
toring) 2.0 m telescope on Haleakala has been carrying out 
observations of A GN and other variable objects (including 
GRB s) since 2001 (|Yoshii et al.ll2002j . l2003l : lKobavashi et al] 
2003). The telescope is equipped with dual optical and in- 
frared channels, allowing simultaneous observations in two 
bands. 

We initiated MAGNUM observations starting at 06:59 
UT, acquiring a sequence of dithered exposures over the next 
~2 hours in a large number of filters, including RI in the 
optical channel and YJHK in the infrared channel. The 
MAGNUM FOV is small, and generally only one star was 
present in the field and away from the chip edge at all dither 
positions. Therefore only a single star was used to establish 
the calibration in each filter. In the H, K, and K-band ob- 
servations the star at q=355. 066002, (5=31.793428 was used 
for this purpose; for R, I, and J the star at q=355. 058815 
(5=31.780569 was used. The second F-band exposure unfor- 
tunately contained no usable reference star. However, com- 
parison of exposures in other filters and at other points in 
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the night suggest that conditions were photometric, and so 
calibration was achieved by comparison to the first Y-band 
exposure (with a small aperture correction.) 



2.8 Kuiper Observations 

Shortly after the GRB trigger, we initiated imaging obser- 
vations at the 1.54m Kuiper telescope, operated by Steward 
Observatory and located on Mt. Bigelow. Observations be- 
gan at 04:37:08 UT and continued until 08:55:41 UT, mostly 
in the R and I filters with some additional observations in V. 
Images were reduced and combined in IRAF using standard 
techniques. The /-band images were not dithered and so an 
archival fringe frame was used to subtract the fringe pat- 
tern. Photometry was performed in IRAF using secondary 
standards. 



2.9 Late-Time Afterglow Observations 

To try to constrain the late-time (t > 12 hr) behavior of 
this burst, additional follow-up was carried out on the 3.6-m 
New Technology Telescope and at GROND. We observed the 
burst location on NTT using the infrared imager SOFI in a 
series of J, H, and /("-band exposures, and additionally in H- 
band only on the following night. Photometry was calibrated 
relative to our IR secondary standards. 

GROND is a seven-channel instrument which has been 
mounted on the ESO 2.2-m telescope at La Silla, Chile, since 
April 2007. GROND began observations of GRB 071025 on 
2007 Oct 26 at 01:50 UT and completed one 8 minute ob- 
serving block and two 20 minute observing blocks. In total, 
9 images were taken in the g'r'i'z' bands and 216 were taken 
in the NIR. Each NIR image was 10 s long; the optical im- 
ages varied in length from 137 to 408 s. The images were 
reduc ed using the GROND pipeline (|Kupcii Yoldas et al.l 
2008), with all images combined into a single stack for each 
filter. For consistency with other measurements, photometry 
was performed using aperture photometry calibrated to our 
secondary standards in JHK. For g'r'i'z' bands, images are 
calibrated directly relative to spectroscopic standard stars 
SA 114-750 and SA 114-656. 

Poor agreement is observed between the NTT and 
GROND observations (and between the overall SED at this 
time and earlier data) using a standard l"aperture, even 
though these epochs are effectively coeval. We have re- 
examined these data and find no clear evidence of problems 
in the reduction or photometry, although the afterglow ap- 
pears extended in the N-S direction in the GROND iif-band 
frame, suggesting that it might be blended with a nearby 
source or image artifact. No neighboring source is observed 
in the NTT imaging, and the deep Keck optical imaging 
shows no object within ~3" of the afterglow position ( H2.10p . 
However, to guard against this possibility we performed the 
photometry in the GROND J and H channels and all NTT 
channels using a small aperture (0.5") in all bands. This 
smaller aperture provides good consistency between the two 
observations and is used in our analysis. 



2.10 Keck Observations 

To help rule out a low-redshift origin for this bu rst we im- 
aged the field around GRB 071025 with LRIS l|Oke et all 
Il995h on the Keck I telescope on 2008 Aug 02 using the 
g and R filters simultaneously under excellent conditions. 
Total exposure times were 1050 s in ii-band and 1140 s 
in <?-band. Consistent with the large photometric redshift 
inferred from the SED, no significant flux due to a host 
galaxy was detected at the location of the optical/IR after- 
glow (the nearest object is a pair of faint point-like sources 
located 3" to the northeast). Forced photometry at the po- 
sition of the optical afterglow, calibrated using unsaturated 
secondary standards, gives a limit (3-cr) of R > 26.5 mag, 
g > 27.0 mag. 



2.11 Field Calibrations 

To improve upon the photometric accuracy of 2MASS, we 
stacked together all observations of the GRB field acquired 
by PAIRITEL during the night of 2007 October 25 UT and 
calibrated a set of isolated, high-S/N stars present in the 
field in all or nearly all dither positions relative to 2MASS. 
These magnitudes were used in place of 2MASS magnitudes 
directly and are presented in Table [S] 

For the optical filters, on the night of 2009 June 19 
we observed the field of GRB 071025 using the Nickel lm 
telescope at Lick observatory. Conditions were photometric 
throughout the night. Three exposures were acquired in R- 
band and one each in I and g band and stars within the field 
were calibrated by com parison to repe ated observations of 
PC 1633 and PG 2336 (|Landoltlll992T ) at varying airmass, 
calibrating reference stars within the field. A second calibra- 
tion was conducted on 2009 Sept 28 using repeated obser- 
vations of standard fields PG 1633, PG 2336, PG 0213, and 
SA 110; the results were found to be completely consistent 
with the June calibration. 

No field calibration was performed in the Y-band. To 
calibrate the observations in this filter, we derived our own 
transformation equation for calculating Y magnitudes of ref- 
erence stars given photometry in nearby bands by fitting a 
simple linear regression model to the photometry available 
online at the UKIRT webpagfl (The Y ~ J colour was fit 
as a linear function of J — H, and the residuals were then fit 
to a linear function of I — J.) The transformation equation 
Y = J + 1.104(J - H) - 0.11(7 - J) - 0.03 was found to 
accurately describe the observed Y-band magnitudes for the 
available standards (with photometry in all four bands) with 
an RMS of < 0.03 mag. We therefore applied this equation 
to calculate the Y magnitudes for secondary standards in the 
GRB 071025 field using the calibrated IJH photometry. 

The final calibrated magnitudes for these stars are pre- 
sented in Table [6] 



5 http://www.jach.hawaii.edu/UKIRT/astronomy/calib/phot_cal/fsjzyjhklm.dat 
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Table 1. Light Curve Best-Fit Parameters 



Parameter 


symbol 


value 


CI rising index 


«l,r 


-1.66 ± 0.15 


CI fading index 


£*l,f 


1.73 ±0.21 


C2 rising index 


«2,r 


-11.0 ± 2.1 


C2 fading index 




1.27 ± 0.04 


CI peak time (s) 


*pk,l 


575 ± 42 


C2 peak time (s) 


*pk,2 


1437 ± 17 


Ratio of C2/C1 peak flux 


F 2 


0.24 ±0.03 


Colour change across CI peak 


A/3 1(rf) 


-0.20 ±0.14 


Colour change between C1/C2 


A/3i2 


-0.26 ± 0.12 


Flux at t = 10000s 


Fr 


5.76 ± 1.15 




Fi 


34.4 ± 3.48 




Fy 


84.45 ± 8.66 




Fj 


118.9 ±4.79 




Fh 


155.4 ±6.26 




Fk 


250.4 ± 15.1 



Summary of free parameters fit in the light curve model. 
Peak times are for the J-band filter. Flux parameters are not 
corrected for Galactic extinction; uncertainties include added 
systcmatics. "CI" refers to the first light-curve component; 
"C2" refers to the second component. 



3 ANALYSIS 

3.1 Early-Time Afterglow Evolution: Rise, Fall, 
and Reddening 

All photometric observations of GRB 071025 during the first 
night are presented in Figure [T] Several features are immedi- 
ately apparent. First, the afterglow was caught during what 
appears to be its initial optical rise, brightening by ~1.5 mag 
from the first detections to the peak in all filters. Second, 
the evolution is not single-peaked: a limited rebrightening 
is observed at ~1800 s. Third, the burst is extremely red, 
with R~K ~ 6.5 mag. Finally, no dramatic colour change is 
evident. This is not to say that there is not finer-scale colour 
evolution, however — as will be discussed later, the best-fit 
curves plotted in Figure [T] correspond to a chromatic model 
which is shown to produce a large improvement in \ 2 rela- 
tive to the monochromatic case. 

The empirical model used to fit this burst is describe d 
in previous work jPerlev et alj l2008bl ; iBloom et ai]|2009l) . 
In brief, our method fits al l data in all filters s imult aneously 
using a series of summed iBeuermann et all (119991 ) broken 
power-law functions, in this case two per filter. The sharp- 
ness parameter s was fixed at 1 (allowing it to vary resulted 
in insignificant improvement to \ 2 ). For each component, 
the pre-break (rising) power-law index a T is constrained to 
be the same for all filters, as is the post-break (fading) in- 
dex Of. The overall flux-normalization factor in each filter 
is arbitrary, determined by the best fit to the data. Colour 
change is modeled by assuming that any overall change takes 
the form of a variation in the intrinsic spectral index f3, i.e., 
F„,2 = F Vl \ x (^■) Af912 . Colour is allowed to vary between 
components (simply using the above equation to tie the nor- 
malised flux parameters of each component) and between 
rising and falling segments of an individual component (by 
allowing the break time of the Beuermann function to vary 
as a function of frequency, using the central frequency of 
each filter). It should be emphasised that this method makes 




Figure 1. Early-time multiband optical and infrared light curves 
of GRB 071025 fit to our empirical light curve model. The after- 
glow is caught during its rise at ~30 s in its rest frame (assuming 
2 = 5), and exhibits a double-peaked structure before fading again 
as a simple power-law. The RAPTOR unfiltered data has been 
shifted to match the I-band data. Magnitudes are Vega-based and 
not corrected for extinction. 

no assumptions about the overall SED, since only the vari- 
ation in P is constrained. Indeed, the fitting method can 
be used to generate a best-fit observed SED in all available 
filters, using all available data, at any chosen time ( i]3.2|l . 

Two components (best-fit parameters are summarized 
in Table [TJ are found to provide an excellent fit to the 
dataQ The light curve brightens quickly between our first 
detections at 180 s with a power-law of approximately 
Qi, r = —1.66 ± 0.15 to a smooth peak at 580 s, then fades 
until about 1200 s. At that point the afterglow briefly re- 
brightens, peaking again at ~1400 s before fading as a sim- 
ple power-law («2,f = 1-27 ± 0.04) for the remainder of our 
observations. The \ 2 residual, assuming no colour change, 
is 222.4 per 154 degrees of freedom (dof). Permitting colour 
change improves the fit significantly: allowing the param- 
eters A/3i2 (describing the change in intrinsic spectral in- 
dex between the first and second component after peak) 
and A/?!^) (describing the change in intrinsic index be- 
tween rising and falling portions of the first component) to 
both vary, x 2 /dof improves to 197.4/152 which (according 
to the /-test) is significant at > 99.9% confidence. Most of 
this change is associated with the transition to the second 
component (A/3i2 =■ —0.26 ± 0.12, versus a not-significant 
colour change across the first peak of A/3 1 ( rf ) = — 0.20±0.14) 
but unfortunately, although the need for overall red-to-blue 



6 The first two Super-LOTIS points are an exception, both of 
which deviate from the fitted model by 2—3 a. Given the low 
signal-to-noise detections and large degree of time-binning in both 
cases, these points are not included in the fit, although the low 
flux observed in the second observational window, which covered 
the peak of the light curve, is nevertheless surprising given the 
behavior in all three PAIRITEL bands and in RAPTOR data at 
that time. 
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Figure 2. The gamma-ray (Swift BAT) and X-ray light curves 
(Swift XRT) of GRB 071025, compared to the J-band light curve 
out to late times. The X-ray light curve is rapidly fading during 
the optical-IR rise, probably due to high-latitude prompt emission 
(the light curve connects smoothly with the BAT light curve at 
these times if scaled to the X-ray flux, as shown.) Both optical 
and X-ray light curves fade with an unbroken decay at late times 
but with different decay slopes: a op t = 1.27 ± 0.04 versus ox = 
1.56 ± 0.03. 

colour change is clear, its nature cannot be clearly distin- 
guished by this methodology. We will further examine sce- 
narios for this possible colour change in H4.2\ 

The X-ray light curve (Figure [2} was fit using a similar 
method (but with only a single "filter" , simplifying the pro- 
cess significantly). Again, two summed functions are found 
to provide an acceptable fit to the data. However, the first 
component is a rapidly-declining, unbroken power-law with 
Qx.init = 3.1 ± 0.2. This initial segment connects smoothly 
with the BAT prompt em ission, as has been s een for a large 
majority of Swift bursts (|0'Brien et alj [20061 ). The optical 
peaks unfortunately fall during an orbital gap in the XRT 
coverage, but by the end of the first observations the power- 
law is already clearly flattening, almost certainly due to the 
transition from the rapid decay phase (O'Brien et af]|2006l ) 
to a standard afterglow dNousek et al.l 20061 ). Coverage re- 
sumes approximately an hour later, by which time the X- 
ray light curve is fading rapidly in an unbroken decay with 
ax = 1.56 ±0.03. 



3.2 SED and Photometric Redshift 

At 10000 s after the burst, the evolution of the light curve 
has given way to a simple power-law decay dominated by 
only a single component. Moreover, thanks to the MAG- 
NUM observations, photometry is available in all colours 
within a relatively short time span surrounding this epoch 
with high S/N in JHK. We therefore choose this time as 
the extraction point for the overall spectral energy distribu- 
tion (SED) of this burst, using the model fluxes from our 
fit as described above. (These fluxes are consistent with the 
MAGNUM and PAIRITET photometry measured at this 



epoch specifically.) All fluxes are corrected for Galactic ex- 
tinct ion (relatively small at Eb-v = 0.07 mag in this direc- 
tion; ISchleeel et aflll99Sh 

The la uncertainty on the fit parameter was combined 
in quadrature with an estimate of the calibration uncertainty 
in each filter. In the J and H filters, where the afterglow is 
comparable in colour to reference stars (which show negligi- 
ble scatter), we use an uncertainty of 0.04 mag; in K where 
the afterglow colour is much redder than our reference stars 
we conservatively increase this to 0.06 mag. This incorpo- 
rates both the absolute and relative calibration accuracy of 
2MASS (estimate d at ~0.02 and 0.011 mag, respectively; 
ICohen et aT1l2003l and 2MASS online documentation^]), ef- 
fects of variation of the effective wavelength A e ft from its ref- 
erence value due to a non-standard spectrum (< 0.02 mag), 
the possibility of strong absorption from ISM or IGM lines 
(very likely < 0.02 mag), and uncertainties in the extinction 
correction (< 0.01 mag). In Y"-band, we use an estimate of 
the photometric scatter of the high-S /N REM reference stars 
to the interpolated secondary standards (0.1 mag). We also 
use 0.1 mag in /-band due to the redness of the afterglow in 
this band and the possibility that Lyman-a may be affecting 
the flux towards the blue filter edge if the redshift is z > 5.0. 
In i?-band a large uncertainty of 0.2 mag is used, although 
because R is almost certainly heavily blanketed by the Ly- 
q forest we use this filter only to place a limiting value on 
the redshift and exclude it from fits to the extinction pro- 
file. The resulting SED (fit with various models, explained 
below) is plotted in Figure [3] 

The sharp dropoff towards the i?-band is suggestive of 
high redshift. However, the spectral slope observed even well 
redward of this apparent break is quite red (/3 ~ 1.64, as 
shown by the dashed straight line in Figure [3J, suggesting 
that significant extinction is likely present as well. In order 
to quantitatively constrain the redshift z, we fit the data set 
with a large number of different extinction models (detailed 
in N3.3P at varying redshifts. Absorption due to the Lyman- 
a forest is taken into account using a simple model of the 
average opaci ty of the IGM as a function of z and A from 
iMadaul (|1995l ). The extinction column Ay and the spectral 
index ft were constrained to be positive: negative extinc- 
tion is unphysical, while a negative spectral index would be 
both much bluer than any previously observed aft erglow and 
in di sagreement with standard afterglow theory (I Sari et al.l 
Il998h . 

The HIRES spectrum discussed in lFvnbo et al.l (|2009T I 
shows a trace extending from the limit of the spectral range 
at 7950 A down to 7550 A, blueward of which no flux is 
detected. While the quality of this spectrum is poor, the 
nondetection of Lyman-a puts a robust upper limit on the 
redshift of z < 5.2, so this was treated as the maximum 
redshift. Regardless of the extinction law, no known dust 
curve is able to reproduce the extremely steep I — R slope 
without invoking Lyman-a blanketing of the i?-band, which 
becomes significant at z ~ 4.0. Even after including a vari- 
ety of extinction templates (below), the lower limit on the 
redshift (95% confidence) is z > 4.4. Treating redshift as a 
free parameter, the best-fit z is dependent on the extinction 
law but is approximately z = 4.8 ± 0.2 (lcr). In the remain- 



7 http://www.ipac.caltech.edu/2mass/releases/allsky/doc/sec4_8.html 
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Figure 3. Spectral energy distribution of GRB 071025 inferred 
from our broadband photometry, fit with different extinction 
models. Note the spectral flattening between J and H that con- 
trasts with red H — K and I — Y colours. (The R— I colour is due 
to absorption by the Lyman-a forest.) Traditional models (such 
as SMC-like extinction, shown here as a dotted line) cannot re- 
produce this feature and give poor fit residuals (x 2 /dof= 20.8/2). 
The supernova-dust model of Maiolino et al. (2004), shown as the 
solid line, is an excellent fit ( x 2 /dof= 0.81/2). The dot-dashed 
line represents the intrinsic afterglow SED (for the SN model) 
without extinction applied but including our model of the IGM 
opacity at this redshift. 



ing discussion we will assume a fiducial value of z — 5.0; 
however, similar conclusions apply to other redshifts within 
the constrained range (4.4 — 5.2). 



3.3 Extinction Profile 

Qualitatively, the SED presented in Figure [3] is unusual 
among GRB afterglows due to the presence of an apparent 
inflection: while the K — H and Y — I colours are very red, 
between H and J bands the slope is quite flat. This flat- 
tening is quite significant (e.g., the U-band point is more 
than 0.2 mag below an interpolation between J and K) and 
suggests that the afterglow of GRB 071025 is subject to 
a complex reddening profile. To try to distinguish differ- 
ent possible models, we therefore fit many different extinc- 
tion laws to the photometric SED, including Milky Way, 
L MC, and SMC curv es estimated using the parametrization 
of iFitzpatrickl {l999) as implemented in the GSFC IDL as- 
tronomy user's libra ry, the starburst-galaxy Calzetti curve 
dCalzetti et al.ll2000[ ), and the high-z QSO extinction law 



is free but limited to be (3 > 0. A summary of the goodness- 
of-fit x 2 for each fit model is presented in Table [2] 

A large family of models, including the Milky Way 
and Large Magellenic Cloud curves as well as the ex- 
tinction curve s derived from a f ew recent highly red- 
dened GRBs jKriihler etHI 120081 : IProchaska et all 120091 : 



from Maiolino et alj (|2004l ). The intrinsic spectral index ft 



lElfasdottir et ah" 20091 ) . display a prominent 2175 A bump. 
We can strongly rule out such a feature: at the observed red- 
shift, the broad absorption signature would fall in or near 
the J-band. Formal fits using these extinction templates (re- 
gardless of Rv) return Ay of zero in all cases (our fits do 
not permit negative extinction.) 

An SMC-like extinction curve (dotted line in Figure [3| 
provides a visually reasonable- looking fit to our data, but the 
X 2 /dof is unacceptable at 20.8/2. This is again no surprise: 
the SMC extinction curve increases rapidly and monoton- 
ically with decreasing wavelength and cannot produce the 
flattening in our SED. The featureless Calzetti law similarly 
produces a poor fit because it cannot produce the deviations 
from a power-law evident in the photometry. 

We also attempted a gener al fit using the full param- 
eterization of IFitzpatrickl (Il999f ). but even if the 7 and xo 
parameters of this model are fixed and the ci an d Ry pa- 
ramet ers are tied C2 using e.g. the correlations of iReichartl 
l|200ll ), the solution is underdetermined. If the intrinsic spec- 
tral slope f3 is fixed, the solution is exactly determined; for 
e.g. p = 0.65, we derive R v = 5.26 ± 0.53, c 2 = 0.17 ± 0.12, 
c 3 < 0.2, c 4 = 1.03 ± 0.32, xV dof = 1-49/0. However, this 
combination of parameters (small C2 and low or zero C3, 
indicating a shallow near-UV extinction law and negligible 
2175A bump) is unlike any sightline in the local universe 
observed to date. We also fit th e data to the general extinc- 
tion curve of lA. Li et al] (|2008l ). fixing /3 = 0.65 and c 4 = 
to avoid underdetermination, but the ci and C3 parameters 
did not converge. 

However, one previously observed extinction law per- 
forms extremely wel l at m atching the observed fea- 
tures. iMaiolino et"all l|2004l ) presented observations of 
the reddened z = 6.2 broad absorption line quasar 
SDSSJ104845. 05+463713, comparing NIR spectroscopy of 
the source to optical spectra of low-redshift quasars of the 
same class to estimate the extinction law. The inferred curve 
of this object is notable for a distinct flattening between 
1800 — 3000A, and was interpreted (and modeled quantita- 
tively) by that paper as the signature of dust synthesised 
in supernova explosions. We fit a polynomial to the solid 
(Z = 10~ 4 Z e ,M = 25M Q ) curve displayed in Figure 2 of 
that paper and used the resulting extinction curve to fit 
our observed photometry[f] The result is an excellent match 
(x 2 /dof = 0.81/2) and is shown as the solid line in Figure 
[3] The associated extinction column is A3000 = 1-09 ± 0.20 
mag0 

The best-fit value of the intrinsic spectral index /3ir. 
as inferred at the SED extraction epoch is /3m = 0.94 ± 
0.14, quite typical of other afterglows at this stage. This 



8 Our i^-band point is not covered by this figure, as the cor- 
responding rest wavelength is shifted out of the IR window at 
z = 6.2. We assume an approximately linear extinction law in 
1/A below Arest < 3300 A. 

9 The Maiolino extinction curve is normalised to A3000 instead 
of Ay (the V-band at 2 > 5 is shifted into the mid-IR). 
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Table 2. Results of Extinction Fits 



Dust Model 


P 


A v 


Ry 






mag 




none 


1.64 ±0.08 







SMC 


0.08 ±0.42 


0.12 ± 20.6 


2.73 


MW 


1.64 ±0.08 


< 0.07 


3.1 


LMC 


1.64 ±0.08 


< 0.07 


3.2 


GRB080607 


1.64 ±0.08 


< 0.12 


4.0 


Calzetti 


0.00 ±0.80 


1.42 ± 0.68 


4.0 


Fitzpatrick 


0.65 


2.52 ± 0.97 


5.26 ±0.53 


Maiolino SN 


0.96 ±0.14 


1.09±0.20 a 





number of calibration stars detected at high S/N, the sta- 
tistical errors on the zeropoint are small. Possible sources 
of systematic uncertainty (beyond the minor effects we have 
already discussed and included) we have considered include: 

Instrumental colour terms. PAIRITEL uses the 
same telescope, filter, and camera system as the 2MASS 
survey, and so there is no reason to expect colour terms as- 
sociated with the optics to be present. However, the presence 
of a significant bandpass difference could cause systematic 
discrepancies in calibration relative to field stars (see also 
i|2.6[) . in particular in H and K bands where the afterglow 
colour is much redder than any of the bright stars used for 
calibration. We inspected the magnitudes derived from stars 
in our deep stack as compared to the stars in 2MASS to 
search for a correlation between the magnitude offset and 
colour; none was found. 

Strong atmospheric variations in the effective fil- 
ter bandpass. The infrared absorption bands associated 
with water in Earth's atmosphere exhibit time-variability, 
even within the observational windows. In particular, the 
exact shape of the J-band transmission f unction depends on 
the amount of precipitable water vapor (I Cohen et al. I [20031 ; 
however, the effect is small, with less than 2% variation in 
relative magnitudes), and the H-hanA contains a water ice 
absorption band which could introduce similar variations. 
Time-dependent absorption may therefore introduce tempo- 
rary colour terms not evident in the all-night stack. There- 
fore, we carefully inspected the time evolution of the ob- 
served zeropoints in all three bands. A small amount (up to 
0.3 mag) of total transmission variability is indeed observed 
during the first 20 minutes, after which the zeropoint in 
all three bands is nearly constant within uncertainties. No 
significant variation is observed in the difference between 
zeropoints in different PAIRITEL bands, nor is any corre- 
lation observed between the overall zeropoint and the dif- 
ference in zeropoints between two bands that would suggest 
chromatic variations in the transmission. Furthermore, the 
zeropoint appears constant (within our uncertainties) after 
~1400 s (the SED is determined at 10000 s). The MMTO 
cloud camera shows no evidence of significant cloud cover 
at any point during the night, and weather archives indicate 
warm and stable conditions during the observation. Further- 
more, in addition to PAIRITEL (Arizona), the Lick J-band 
(California) and MAGNUM (Hawaii) coeval measurements 
both give consistent results for the infrared magnitudes, giv- 
ing additional confidence in our results; in particular both 
PAIRITEL and MAGNUM JHK data sets show the puta- 
tive extinction feature independently. Therefore, we have no 
reason to believe that our SED is significantly affected by 
absorption features in Earth's atmosphere. 

Intrinsic deviation of the GRB spectrum from a 
power-law. We have assumed in our fits that the intrinsic 
spectrum of the GRB was a simple power-law, as generally 
predicted by synchrotron theory. This assumption could, in 
principle, be violated. However, the most natural deviation 
from a power-law SED that might be expected (a spectral 
break within the optical/IR band) would create downward 
curvature in the intrinsic SED and actually require addi- 
tional dust to produce the upward inflection feature that is 



Summary of key parameters from fits of various dust models 
to the SED of GRB 071025 as modeled at t=10000 s. A 
redshift of z = 5 is assumed in all cases. 

a Value is ^3000- 

value is also consistent (albeit only marginally) with the 
theoretically expected value based on the observed X-ray 
spectral index (intrinsic /?x = 1.15 ±0.12) if a cooling break 
is present between IR and X-ray bands: in this case /3ir = 
/3x — 0.5 = 0.65 ± 0.12.) Imposing this constraint as a prior 
on the fit to Av and /3m,, we measure Ay = 1.27±0.20 mag. 

Alternatively, the SED is also consistent with the pres- 
ence of no cooling break: at the extraction epoch the com- 
bined IR-through-X-ray SED is well-fit (x 2 / dof = 1-14/3) by 
a single power-law with /3ir,x = 0.88 and Av = 1.19 ± 0.20 
of Maiolino dust, both consistent with the values inferred 
from the optical data alone. However, the X-ray flux at this 
time is clearly fading faster than the optical light curve (Fig- 
ure!]]): if this is not due to the presence of a moving spectral 
break such as a cooling break, the spectral index itself would 
have to be slowly evolving (implying evolution in the elec- 
tron index p) . 

In support of our general conclusion of a significant 
amout of dust extinction, we note that a large amount of 
absorption is inferred from the X-ray spectrum also: we mea- 
sure an equivalent column of Ah = (3.2 ± 0.8) x 10 22 cm~ 2 . 
Although the scatter in the ratio of Av / Nh for Swift bursts 
is nearly an order of magnitude, using the average value from 
ISchadv et all \200H ) this column corresponds to an extinc- 
tion of Ay ~ 4 mag. 

3.4 Further Investigations of the IR Calibration 

The inferrence of SN-type dust for this object depends sen- 
sitively on the accuracy of our photometric calibration, and 
statements of its significance relative to the SMC fit depend 
equally critically on the precision in the JHK bands being as 
good as we claim: the Maiolino model is no longer preferred 
at > 95% confidence if, for example, additional uncertainty 
of >0.075 mag (in addition to the systematic uncertainties 
already applied; ij3.2[) is added in quadrature to all SED 
data points, or if >0.1 mag is added to just the ff-band 
point (dependence on the other data points is much more 
robust: an addition of >0.2 mag to K is required, and any 
one of the J, Y, or I points could be removed completely). 
Therefore we have scrutinised in detail our infrared calibra- 
tion procedures with particular emphasis on the PAIRITEL 
data. Because of the large number of exposures and large 
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Figure 4. Time-dependent spectral energy distribution of GRB 
071025 inferred after dividing the data into six different windows 
and re-fitting the flux parameters at each epoch using the light 
curve model. The resulting SED is then fit for spectral index 
and extinction column A3000 a t each epoch individually using a 
Maiolino extinction profile. The characteristic flattening between 
J- and //-bands is observed at every epoch (except at late times, 
when photometric errors are large) with no significant variation in 
its strength, increasing our confidence that it is a feature extrinsic 
to the GRB. 



observed. An SED modeled as the sum of two components 
(a steep power-law dominating K— band and a shallower 
power-law dominating J-band) would produce upward cur- 
vature, but cannot reproduce the sharpness of the observed 
feature unless the spectral index of the steep power-law is 
unrealistically red (/? > 4). Additionally, it would be surpris- 
ing that both components would rise and fall in synch with 
each other throughout the complex early evolution of the 
light curve, as is observed. Indeed, evidence for Maiolino- 
like dust is observed at every epoch with no significant vari- 
ation in its strength or wavelength (Figure 3] also >i4.2|l with 
the exception of the final (GROND+NTT) SED, when the 
photometic uncertainties are too large to place any strong 
constraints on the extinction law. 

Absorption from a DLA host system or Lyman- 
q forest. If the host galaxy is at the maximum redshift of 
.z=5.2, our mean-opacity model of the Lyman-a forest may 
significantly underestimate the impact of hydrogen absorp- 
tion on the /-band. To represent the most extreme possible 
case, we reran our dust models after adjusting the /-band 
flux upward by 20% (the approximate maximum dimunition 
expected in the Kuiper /-band filter assuming 100% opac- 
ity blueward of 7550 A, the limit on any DLA imposed by 
the HIRES spectrum) at z = 5.2. Even in this case, the 



Maiolino dust profile is strongly preferred (x 2 /dof = 2.5/2, 
versus 11.9/2 for SMC-like dust.) 



4 DISCUSSION 

4.1 Rise of the Forward Shock and Constraints on 
the Lorentz Factor 

The nearly achromatic first peak in the light curve shows the 
major hallmarks of the initial rise of the afterglow due to hy- 
drodynamic deceleration of the fireball: a steep rise with no 
significant evidence of colour change across the peak. Al- 
ternative possibilities can be generally ruled out: for exam- 
ple, the transition of the synchrotron peak frequency (which 
would also produce a peak were it to occur after the ini- 
tial deceleration) would rise slowly and exhibit a blue-to- 
red colour shift of A/3 = (p - l)/2 - (-1/3) = p/2 - 1/6, 
completely incompatible with the observations of no colour 
change or even limited red-to-blue evolution at this time. If 
the peak were due to dust destruction we would also expect 
significant colour change during the rise itself, which is not 
apparent in the data. (We will examine the possibility of 
dust destruction in more detail in i]4.2[) . 

Within the category of hydrodynamical effects, there 
are then three possibilities for the rise of the afterglow: peak 
of the reverse shock, peak of the forward shock, or an off-axis 
jet. 

We will consider the jet model |Granot et al.l 120021 : 
lGranotll2005l ) first. In this case, the outflow is assumed to 
be strongly collimated with an observer located outside both 
the jet opening angle 9 (observing angle o bs) and Lorentz 
cone 1/F (for a uniform jet; the theory can be suitably m od- 
ified for a structured outflow: iKumar fc Granotl 120031 ) . As 
the jet decelerates, a peak in the light curve will be ob- 
served once the flow has decelerated sufficiently for the 1/F 
cone to expand past the observer line of sight; this model has 
shown reasonable su ccess representing the rising light curves 
of e.g. XRF 080330 (iGuidorzi et alj|2009h and GRB 080710 
(Kriih ler et al.ll2009bl ). However, we are disinclined to favor 
this model on the grounds that it is expected to produce a 
very rapid post-break decay (a > 2), which is not observed 
at late times (ct2,f = 1-27 ± 0.04). This could be accounted 
for by associating the second component (which dominates 
the late-time deca y) with an on-axis wid e jet undergoing its 
initial rise (as in jKruhler et alj|2009bh . but this model is 
somewhat contrived in our case, requiring fine-tuning of the 
physical properties of the two jets to accommodate the large 
variation in their jetting times while still ensuring that they 
peak within a factor of ~2 in time and flux. Alternatively, 
refreshed shocks and continuous energy injection out to late 
times could also be invoked to explain the two-peaked struc- 
ture and lack of late decay within this model. Even in that 
case, another criticism of this model is that the isotropic en- 
ergy release observed for this burst (Ei so = 6.5 x 10 53 erg) 
is not expected for a burst seen off-axis. 

Next, we c onsider if the initia l rise could be due to the 
reverse shock (|Sari fc Piran| [l999). This model is particu- 
larly attractive, as the overall light curve qualit atively looks 
impre ssively similar to the theoretical curve of IZhang et al.l 
(2003): the first peak corresponds to the reverse shock and 
the second peak to the forward shock. However, the initial 
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rise is somewhat slower than expected from simple analytic 
models. The assumed reverse shock rising index a\, r de- 
pends on the assumed zero time to (which was set to the 
trigger time in the above fits), but to would need to be 
shifted back in time by an amount greatly in excess of the 
duration of the burst itself to match the predicted f 3p-5 / 2 
predicted for the reverse shock rise in the slow-cooling case 
l|Kobavashill2000h . The alternate fast-cooling case predicts 
a slower rise (too slow: t 13//16 ) and also a bluer spectrum 
than is preferred by our extinction modeling. A wind model 
also requires fast-cooling and a blue spectrum, and an even 
slower rise (t 1/2 ). Therefore a reverse shock is not our pre- 
ferred paradigm either, though we are hesitant to rule it 
out on the basis that the known complexity of early after- 
glows and the failure of even late-time closure relations to 
properly predict the decay rate a (e.g.. iRvkoff et al]|2009l ) 
suggest that the quantitative details of light curve behavior 
may not be an especially reliable way to evaluate different 
models. 

The most straightforward scenario for the initial rise is 
the formation of the forward shock as the burst ejecta decel- 
erate s into the surrounding medium (e.g.. lRees fc Meszarosl 
Il992h . In this case a = —2 for v < v c , which is still somewhat 
too fast but still consistent with the data within 2a if to is 
moved backwards in time by about 30 s. In this model, the 
second peak is presumably due to additional energy input 
from the central engine into the forward shock, perhaps in 
the form of a slow-moving shell that catches up at around 1 
ks (|Rees fc Meszarosl 1 19981 ). This model is generally consis- 
tent with all available observations including the apparent 
rapid rise of the second component, though the observed 
significant (albeit minor) colour change is not predicted. It 
could be due to the passage of a cooling break (though would 
imply v > v c initially and a too-steep a = —3 during the 
rise) or another effect such as variation in the electron index 
P- 

Interpreting this feature as a forward shock enables us 
to meas ure the ini t ial Lo rentz factor of the e x plosio n. Follow- 
ing e.g. iMeszarosI (|2006l ) and IRvkoff et al l (|2009h . this can 
be estimated from observable parameters via the following 
relationship: 



2r d 



2( 



560- 



3i5iso 



327mm p cVpk, z ' 



a/8 



3-Biso,! 



7?o.2n ^ k , z , 10 

Here i?iso,52 is the isotropic-equivalent energy release in 
units of 10 52 erg, 770.2 is the radiative efficiency in units of 0.2, 
no is the circumburst density in units of cm -3 , and tpk.z.10 
is the afterglow peak time as observed at the burst redshift 
z in units of 10 s. For GRB 071025, using E iso = 6.5 x 10 53 
erg from our spectral model of the BAT data (at z = 5), we 
derive: 

To ~ 20677 "2 /8 n" 1/8 

Compared to direct pair-opacity lower limits inferred 
by the Fermi LAT ijAbdo et al.ll2009al rbh. this is a relatively 
low value of V. However , it is fairly typical of afterglow- 
inferred values (100-1000. IRvkoff et al. 2009 UMolinar7 et al.l 
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Figure 5. Models of colour evolution in the afterglow of GRB 
071025. (a) Infrared and X-ray light curves of GRB 071025 
from PAIRITEL and the Swift XRT showing the divergent be- 
havior in the two bands at these times. The early X-ray light 
curve is probably dominated by prompt emission, which is con- 
tinuing in the BAT as well during this early decay phase, (b) 
The infrared/optical spectral index /?, as measured by a fit to 
PAIRITEL JHK and RAPTOR unfiltered data. Fixed extinc- 
tion A3000 = 1-1 mag is assumed. The SED is observed to redden 
significantly during the obervations. Grey points indicate fits to 
PAIRITEL JHK photometry only, (c) The time-dependent ex- 
tinction column A3000 as measured by PAIRITEL and RAPTOR. 
The spectral index f} is also free to vary in these fits. No evidence 
for variation in the extinction column is observed, ruling out dust 
destruction after ~150 s. 



120071 . lOates et all [20091 . iKriihler et al.ll2009al lbh. This may 
indicate a difference in the types of populations probed by 
the two methods: the intrinsic delay in optical follow-up can 
measure F only for bursts for which the peak is quite late 
(low F), while high-energy photons themselves escape only 
if r is large. Hopefully, in the near future a joint Swift- 
Fermi burst with a luminous afterglow will allow both meth- 
ods for estimation of the Lorentz factor to be compared. 



4.2 Colour Evolution: Limits on Dust Destruction 

Because of the need for vm < v c to explain the slow rise, 
there is no explanation within the standard assumptions 
of afterglow theory for the colour change observed during 
the afterglow. One possible solution would be to invoke a 
time-variable electron index p at early times; a softening 
of the electron distribution during the complex early evolu- 
tion would cause a corresponding softening of the afterglow 
emission. 

Another intriguing possibility, however, is the pho- 
todestruction of dust a long the GRB lin e of sight 
l|Waxman fc Drainel I2OO0I : iDraine fc Haol I2002T ). While we 
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have ruled out this model as being the predominant origin 
of the rise of the light curve based on the modest or ab- 
sent colour change during the rising phase, it is still possible 
that it is occurring on a more subtle level. Because our light 
curve model assumes any colour change is associated with a 
temporal break, it is not clear that such a change would be 
manifest in those models. As a result, we have scrutinised 
the overall colour evolution of this GRB in significant ad- 
ditional detail to search for time-evolution in the extinction 
column Ay. 

The large flat (grey) component of the Maiolino SN- 
type extinction law has the useful feature that the observed 
spectral slope of an SED measured over this region will 
closely match the intrinsic spectral slope even for a large 
extinction column, breaking the degeneracy between the in- 
trinsic spectral index /3 and amount of reddening ^3000 • At 
z ~ 5, the J — H colour (where the extinction law is grey) 
is affected only by the intrinsic spectral index and is nearly 
independent of ^3000, while H — K and J — I are affected 
by both the intrinsic index and reddening. This allows us to 
fit for P and A3000 independently with reasonable reliability, 
even with only a small number of points in the SED. 

We have, therefore, undertaken time-variable extinc- 
tion fits using simultaneous measurements by mosaicing the 
PAIRITEL JHK data to temporally match the early-time 
RAPTOR points, which are a good approximation of the I- 
band (after a small adjustment of —0.08 mag: see i]2.4l and 
Figure [TJ. Dust models were fit to this four-point SED as in 
section i]3.3l Results are plotted in Figure [5] 

As in the case of the complete data set, a Maiolino 
dust model is significantly preferred, with no evidence of 
evolution. In particular, the first mosaic (the only one con- 
temporaneous with bright X-ray prompt emission, which 
is probably the dom inant contributor to dust destruction: 
iFruchter et~ai1l200ll ) gives a modest value of ^3000 = 1.03 ± 
0.31, fully consistent with our measurement at 10000 s. The 
corresponding 95% confidence limit on the decrease in the 
extinction column is A^3ooo < 0.54 mag. As we observe the 
//-band (A e a ~ 2770A in the host frame) rising by at least 
1.5 mag between our first REM and PAIRITEL observations 
and the peak, this clearly rules out dust destruction as the 
cause of the early peak, consistent with our conclusions of 
the chromatic light curve modeling in ij4.ll The entirety of 
the colour variation appears to be due to variation in the 
intrinsic spectrum. 

The significant dust column, combined with the lack of 
variability even during the end of the prompt phase, places 
a limit on the proximity of t his d ust to the GRB. The sim- 
ulations of e.g. iPerna et a l. (2003) suggest that for a bright 
GRB virtually all dust within about 10 pc of the GRB will be 
destroyed, and significant destruction will be observed even 
out to 100 pc. While the exact constraints for this event will 
likely depend on detailed modeling of GRB 071025 specifi- 
cally, this gives an approximate limit on the distance of the 
inferred absorbing dust column from the progenitor of at 
least > 10 - 100 pc. 



5 CONCLUSIONS 

GRB 071025 joins a growing list of gamma-ray bursts caught 
early enough in their evolution to observe the rise and peak 



of the optical afterglow. Interpreting this as the initial rise 
of the forward shock, we estimate T ~ 200 for typical ISM 
densities. The mild red-to-blue colour evolution of the af- 
terglow appears to be due to unknown intrinsic properties 
of the forward shock, rather than dust destruction due to 
irradiation of the burst environment. All of these properties 
are similar to those inferred from early-time observations of 
other GRB afterglows. 

However, the extinction law we measure is nearly 
unique. Most aft erglows with well-ch aracterised SEDs show 
little extinction (|Kann et al.ll2007bT) . and events for which 
significant extinction has been observed have most com- 
monly shown simple SM C-like profiles (e.g., iKann et al.l 
120061 : ISchadv et alj|2007h . characterised by significant cur- 
vature (strong wavelength dependence) but no spectral fea- 
tures. More rarely, featureless or even grey light curves with 
no si gnificant curvature have been inferred for some bursts 
.( e -g-; 



2006; 



Savaglio fc ryiill2004l: Istratta et al.ll2005l : IChen et all 
Perlev et al.ll2008al : |Y. Li et al.ll2008f ). and recently a 



small number of events have been discovered containing the 
clear signature of the 2 175- A bump present in the Milky 
Way and LMC JKriihler et alJliool iProchaska et al.ll2009l ; 



lElfasdottir et al 



2009), though the details of these extinc- 



tion curves show some differences from the average Milky 
Way ISM law. But to our knowledge, no other GRB sight- 
line has shown clear evidence of dust not well-fit either by a 
local extinction template or by a simple, featureless law. 
A possible exceptio n is z=6.3 G R B 05 0904. For this 

(120071 ') favored the 



GRB, an analysis b y IStratta et al 



supernova-type dust of lMaiolino et al 



over standard 

SMC, Milky Way, and Calzetti models. Taken together, 
these two bursts would represent compelling evidence of an 
association between the observed dust model and the chem- 
ical evolution of the universe itself: to date, these events 
are the only bursts at z > 4.5 showing evidence for signifi- 
cant extinction (all other bursts for which useful constraints 
on the extinction law have been possible are at z < 4: 
IKann et al.ir2007bl )Fl 

Strong chemical evolution of the dusty ISM is to be ex- 
pected at z ~ 5 — 6: while most dust at low-to-moderate 
redshifts is thought to have been produced in AGB stars, 
during the first ~1 Gyr following the Big Bang there had 
not yet been time for thes e stars to form in large numbers 
j|Morgan fc Edmunddl2003h . The cosmic age of 1.1 — 1.4 Gyr 
(assuming standard cosmological parameters with ACDM) 
allowed by our photometric redshift suggests that SN-like 
dus10 is still the predominant source of obscuration in 

11 Unfortunately, the case of GRB 050904 is still ambiguous. 
Numerous ot her papers have investigated the du st properties 
of this event jKann et alj l2007al : iGou et all 120071 ') and none of 
these other autho r s pre sented evidence favoring the Maiolino 
curve. Liang & Lj J2009T) have claimed detection of the 2175 A 
feature. Therefore, we downloaded the available data on this 
sourc e llHaislip et alj [20061 : iTagliaferri et al.1 120051: iKawai et all 
2006: lBoer' eTaLll^OOg) and attempted to model the dust profile 
of this event using the same tools applied to GRB 071025, and 
found no evidence for a featured extinction curve. Indeed, the data 
are fully consistent with no extinction at all: our extinction fits 
converged to a simple power-l aw with ~ 1.0, in agre ement with 
the co mprehensive analyses of lKann et al.1 j2007aT ) and lGou et al.1 
ll2007ll . 

12 Although the extinction profile observed is an excellent match 
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galaxies at this epoch and could provide important con- 
straints on the evolution of the firs t galaxies and the pro- 
duction of early dust grains (e.g., iDwek et alJ l2007t >. We 
note that the extinction column inferred from this galaxy 
is even larger than that inferred from the z = 6.2 QSO 
043000 = 0.4 - 0.8 mag, iMaiolino et~ai1l2004h . suggesting 
that even at this epoch, significant amounts of dust are 
present near sites of active star formation. Alternatively, the 
unusual dust could be associated with the relatively nearby 
environment of the GRB only and not necessarily repre- 
sentative of the galaxy itself. However, the survival of this 
late dust limits the distance from the progenitor to at least 
10 — 100 pc, suggesting that its association with the progen- 
itor star-forming region cannot be too close. 

The case of GRB 071025 is illustrative of the potential 
for early-time broadband photometry of GRBs to reveal 
the c hemical history of the early universe (|Hartmann et all 
2009). Well-characterised high redshift bursts are unfortu- 
nately rare (five years into the Swift mission, only five other 
bursts to date have been confirmed to be at 2 > 5), and 
high-redshift events showing significant dust columns are 
even rarer (with the exception of the controversial 050904, 
above, none of the other z > 5 events show evidence for sig- 
nifica n t reddening: Greiner et al.ll2009a ; iRuiz-Velasco et "all 
120071 ; iTanvir et al.1 120091 ; iKann et all |2007bl.) How ever, 
future GRB missions such as EXIST |Grindlavi r2009) arc 
likely to produce a large advancement in our understanding 
of these events. While designed to search for GRBs at the 
redshift extremes (z > 7) and characterise these events 
spectroscopically, infrared photometry and spectroscopy 
acquired of the much more frequent moderate-redshift 
events (z = 4 — 7) will place important constraints on 
the abundance and composition of dust during these early 
stages of cosmic evolution, when galaxies were in the active 
phase of assembly and the first generations of stars led to a 
rapid build up of the metal content of the universe. 
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to the models of d ust produced in supernovae provided by 
IMaiolino et al. we note that this does not demonstrate 

conclusively that this dust was formed in the supernova explo- 
sion itself. Alternatively, the dust could be formed in the ISM 
toraindl2009l ) fro m refractory eleme nts produced in SNe, from 
early carbon stars l lSloan et al.l l2009V or from an unknown pre- 
AGB formation mechanism. 
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Table 3. PAIRITEL Observations of GRB 071025 



^start ^cnd ^exp J H K s 

s s s mag mag mag 



162 


3 


186. 


.9 


23. 


.4 


16. 


160 ±0 


.452 


15, 


933 ± 0. 


.768 


14. 


.098 ± 0. 


413 


197. 


6 


222 


.1 


23 


,4 


15 


448 ± 


.153 


14. 


898 ± 0. 


.198 


13. 


.751 ± 0. 


195 


233. 


.8 


258 


.2 


23. 


.4 


15. 


.486 ± 0. 


.160 


14, 


589 ± 0. 


143 


13. 


397 ± 0. 


138 


269. 


.9 


294 


.4 


23. 


,4 


15. 


300 ± 


141 


14, 


299 ± 0. 


103 


13. 


259 ± 0. 


119 


306 


1 


330 


6 


23 


,4 


15. 


.083 ± 


143 


14, 


.072 ± 0. 


.124 


13. 


.070 ± 0. 


131 


342 


4 


366 


.9 


23. 


,4 


14. 


.814 ± 


126 


13 


.987 ± 


.112 


12 


.905 ± 0. 


136 


378 


6 


403 


.1 


23. 


.4 


14. 


.797 ± 


.117 


13 


956 ± 


.110 


12 


.884 ± 0. 


110 


414. 


.8 


439 


.3 


23. 


.4 


15 


058 ± 


.129 


14. 


.065 ± 


.117 


12 


.854 ± 0. 


112 


451 


1 


475 


.5 


23 


.4 


14. 


624 ±0 


.091 


13 


817 ±0 


.086 


12 


.835 ± 0. 


.098 


487. 


,2 


511 


6 


23 


,4 


14. 


.669 ± 


.084 


13 


.885 ± 


.082 


12 


.606 ± 0. 


.079 


523. 


3 


584 


.0 


46. 


.8 


14. 


600 ± 


.064 


13 


.774 ± 


.064 


12 


.630 ± 


.064 


595 


.7 


656 


.4 


46. 


.8 


14. 


.801 ± 


.083 


13. 


.884 ± 


.071 


12 


.919 ±0 


.080 


668. 


1 


728. 


.8 


46. 


.8 


14. 


.758 ±0. 


.073 


13 


.913 ±0 


.069 


12 


.858 ±0. 


.073 


740. 


.5 


801 


.2 


46 


.8 


14. 


677 ± 


,059 


13. 


871 ± 


.057 


13 


.007 ± 0. 


.073 


813 


.9 


874 


.6 


46. 


.8 


14. 


.799 ± 


.088 


13 


.977 ± 


.083 


12 


.875 ± 


095 


886. 


3 


946 


.9 


46. 


.8 


14. 


.942 ± 


,083 


14, 


.250 ± 


.083 


13 


.136 ±0 


.090 


958 


6 


1019 


.4 


46 


.8 


14. 


982 ±0 


,083 


14. 


.100 ±0 


.076 


13 


111 ±0. 


.095 


1031. 


.1 


1128. 





70. 


.2 


15 


359 ±0. 


,099 


14, 


.466 ± 


.087 


13. 


.417 ±0. 


.094 


1139. 


8 


1236. 


.7 


70 


2 


15 


.374 ±0. 


,078 


14, 


355 ± 


.061 


13 


.355 ±0, 


.075 


1249. 


4 


1346. 


.4 


70. 


2 


15 


.161 ±0. 


,080 


14, 


.384 ± 


.078 


13 


.330 ± 


090 


1358. 


1 


1456 





70. 


.2 


15. 


.247 ±0. 


,066 


14, 


.443 ± 


.063 


13 


,311 ±0, 


.075 


1467. 


7 


1564 


6 


70 


,2 


15. 


.086 ± 


,052 


14, 


.257 ± 


.056 


13 


,296 ± 0. 


.068 


1576. 


3 


1673. 





70. 


.2 


15. 


064 ± 


,052 


14 


.351 ±0 


.061 


13 


,465 ± 


.075 


1684. 


7 


1781. 


6 


70. 


2 


15. 


.176 ± 


,052 


14. 


.405 ± 


.056 


13. 


.435 ± 


.073 


1793. 


3 


1890. 


.2 


70 


2 


15. 


.193 ±0. 


,052 


14, 


.452 ± 


.054 


13 


.528 ±0. 


.073 


1901. 


9 


1998 


.8 


70. 


.2 


15. 


.227 ±0. 


,054 


14, 


.527 ±0 


.059 


13. 


.554 ±0. 


.082 


2010. 


5 


2107. 


3 


70. 


.2 


15 


302 ± 


,063 


14, 


.564 ± 


.066 


13 


.525 ±0, 


.082 


2120. 


1 


2216. 


9 


70. 


.2 


15 


434 ± 


,066 


14, 


.573 ± 


.063 


13 


.739 ±0, 


.094 


2228. 


6 


2325. 


6 


70. 


2 


15 


543 ± 


,078 


14, 


.702 ± 


.066 


13. 


,775 ±0, 


.097 


2337. 


2 


2434. 


2 


70, 


.2 


15. 


637 ±0. 


,066 


14, 


.763 ± 


.066 


13 


,713 ±0. 


.087 


2445. 


9 


2542 


.8 


70. 


.2 


15. 


.644 ± 


,071 


14. 


.927 ± 


.078 


13. 


,887 ±0. 


.101 


2555. 


5 


2652 


3 


70. 


.2 


15. 


.602 ± 


,082 


14, 


.959 ± 


.087 


14 


Ill ±0, 


.125 


2664. 





2760. 


.9 


70. 


2 


15. 


848 ± 


,087 


15. 


.039 ± 


.087 


14 


.006 ± 


.118 


2772. 


6 


2941 


9 


117 





15. 


.985 ±0. 


,094 


15 


100 ±0 


.084 


14 


,081 ±0. 


.116 


2953. 


6 


3196 


.4 


163. 


.8 


16. 


004 ± 


,070 


15. 


.211 ±0 


.072 


14 


,078 ± 0. 


.089 


3208. 


1 


3486. 





187. 


2 


16 


168 ±0 


,075 


15 


.455 ± 


.084 


14 


.405 ± 0. 


111 


3498. 


7 


3812. 


9 


210. 


6 


16 


.314 ±0. 


,082 


15 


420 ± 


.075 


14 


.415 ±0, 


.108 


9108. 


8 


12432. 


.7 


2152 


.8 


17. 


.680 ± 


154 


17 


.196 ±0 


.226 


16 


.074 ± 0, 


.228 


13132. 





16637. 


.0 


2269. 


.8 


18. 


.344 ± 


,228 


17 


.462 ± 0. 


.245 


16. 


,564 ± 0. 


.324 



Time values are measured from the Swift trigger (UT 2007 
Oct 25 04:08:54). Magnitudes arc in the 2MASS (Vega) 
system and not corrected for Galactic extinction. 
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Table 4. Additional Photometry of GRB 071025 



Telescope 


^mid 
S 


Filter 


^exp 
S 


Mag. 






Flux 
fijy 


RAPTOR 


119.5 


clear 


45.0 


> 16.94 






< 608.9 


R APTDR 


9Q0 4 


clear 


900 

zuu.u 


171 07 j_ 

1 i . IO I ZC 


u. 


1 

■ ioy 


/1QPC Q _|_ fi.fi O 

'loo.y in uu.z 


R APTDR 


526.3 




180.0 


16.793 zc 





118 


«QQ PC _|_ 71 Q 

uyo.o zc 1 x.y 


T> APTOR 


739.6 


clear 


180.0 


16.792 ± 





111 


(SQQ 9 4- fi7 Q 
uyy.Zr zc ui .y 


R APTOR 

X V, A X X V_/ X V 


953.3 


clear 


180.0 


16.761 ± 





.112 


71 q 4 _i_ 70 C 

1 X O . '-i _1_ 1 KJ .fj 


R APTOR 


1167.1 


clear 


180.0 


17.517 zL 


Q 


236 


QPCQ ft _|_ 70 1 
OOO. U ZC IU.1 


R APTOR 


1381.1 


clear 


180.0 


17.012 ± 





140 


PC7H Q _|_ ftQ 1 

j 1 u.y zc uy . x 


R APTOR 

X V, A X X vy X V 


1594.9 


clear 


180.0 


17.096 ± 


q 


156 


pc9S 4 4- 70 7 


R APTOR 


1808.7 


clear 


180.0 


17.186 zL 





170 


40ft 4 _|_ 70 c; 

IOU.1 ZC iU.O 


R APTOR 


9099 3 


clear 


180.0 


17.216 zL 





172 


47Q 1 _|_ ftQ Q 

110.1 zc uy .0 


R APTOR 


ooqc; ft 
zzoo.u 


clear 


1 so 


1 7 ftQI _L 

1 i .uyi zc 


ri 

u. 


281 


ouo.o zc uy. 1 


R APTDR 
j.\,ji--r x v/i i 


2448.6 


clear 


180.0 


17.429 zlz 


Q 


227 


QQQ Q _|_ 70 4 

000.0 ZC t O.I 


R APTOR 

Xwi-X X vy IV 


2662.1 


clear 


180.0 


1 o^O 4- 

IO.UUU ZC 


Q 


372 


91 Q C _|_ ftQ 7 

zxy.o zc uo. f 


R APTOR 


2875.7 


clear 


180.0 


18.097 zL 


Q 


413 


91 9 4- fifi ^ 
ZXU.Zr ZC uu.o 


R APTOR 


3124.9 


clear 


240.0 


17.793 zL 


Q 


265 


970 1 _|_ fto 9 

Z. 1 O. X ZC UU. Zr 


Snnor T OTTS 

O UpLl JjU 1 lu 


134.5 


R 


50.0 


> 19.46 






< 59.89 


O Lip CI -JjU llO 


244 2 


T5 
I V 


1 00 


in lon-L 
iy . .Lou zc 


U. 


240 x 


77 ^1 4- 1 ^ ^7 X 
/ 1 .Jl zc 10.0 / 


Snnpr TOTTS 


478.7 


R 


300.0 


19.700 zL 





240 x 


40 09 4- Q c 19 x 


Snnpr T OTTS 


813.7 




300.0 


18.910 zL 





120 


no 40 _|_ 1 40 


Snnor-LOTTS 

k_> Ll 1} \_ 1 Jjv/ X XO 


1315.1 


R 


600.0 


19.520 ± 


q 


.160 


Kft «7 _i_ 7 77 
•j \) • \j 1 _i_ 1 . 1 1 


Snnpr T.OTTS 


1983.3 


R 


600.0 


19.390 zL 


Q 


180 


ftQ OQ _|_ Q 7ft 

uo.oo zc y. 1 u 


REM 


470.0 


Y 


81.0 


lO.U^U ZC 


Q 


240 


1 9PC7 9 _L 94Q Q 


REM 


1281.0 


Y 


181.0 


15.700 ± 


q 


190 


1 1 (57 q 4- 1 qv c: 


REM 


ZrUOO . U 


Y 


332.0 


1 fi 980 4- 

lU.ZrOU ZC 


Q 


330 


ftQ4 ft _|_ 1 70 4 
uoi.u zc x 1 y.i 


REM 


377.0 


j 


91.0 




Q 


200 


1 ^74 9 4- 9^1 9 

IO 1 ZIZ ZOX.Zi 


REM 


1085.0 


j 


181.0 


15.570 zL 





.160 


QQPC C _|_ 1 Oft 4 

yyo.o zc iou.i 


REM 


2304.0 


j 


331.0 


1 cc Qtcn 4- 





110 


19101 4.117c 
xzi x y . x ziz 111 .0 


REM 


185.0 


jj 


82.0 


16.131 ± 


Q 


690 


0,70, 7 J. 17c: 
1 0. t zc 110.0 


REM 


666.0 


jj 


181.0 


1 3 KQO 4- 


Q 


050 


9QA4 3 4- 1 ^9 ^ 


REM 


1623.0 


H 


331.0 


14.380 ± 





080 


1 274 q 4- 1 5Q 


REM 


2795.0 


jj 


81.0 


> 14.90 






< 1161.4 


REM 


97^ n 


K 


89 
OZ.U 


1 3 4^0 4- 


U. 


1 ^0 

. IOU 


9&QQ 9 4- 397 9 

zoyy.z zc ozi .z 


REM 


&qo n 
oyu.u 


K 


212 


1 9 QOO 4- 


n 
u 


070 


ZL793 4- 9QA Q 


Lick 


971 A n 


j 


91 

ZrlU.U 


1 PC Q-l 7 I 

10.01 i zc 




0^0 


709 Q _|_ 91 ft 

/ yz.y zc zi .u 


Lick 


9QQ7 

zyy / .u 


7 

J 


91 
Z1U.U 


1 ^ 0^8 4- 


n 

u. 


0^0 

uou 


70Q Q J_ 1 Q Q 

/ uy .0 zc iy .0 


Lick 


Q97Q n 

oz * y .u 


T 

J 


91 

ZrlU.U 


1ft 1 QC. _L 
-1U. -IOO zc 


u. 


0^0 

uou 


PCQ1 ft J_ 1 ft 1 

oyi.u zc iu. 1 


Lick 


3^9 
ouuz . u 


J 


210.0 


16.227 z\z 





030 


PC4Q ft 4. 1 4 Q 
OIO. u zc 11.0 


Lick 


3846.0 


J 


210.0 


16.378 ± 





030 


47Q 4- 1 9 Q 


Lick 


4129.0 


J 


210.0 


16.504 ± 





030 


491 oj- 11 p. 

IZiX.Zi ZC ll.U 


Lick 


4413.0 


J 


210.0 


1 a KPCQ _L 

1U . iJ ' J O _l_ 


q 


030 


4nn 7 4- 10 Q 


Lick 


iuyo . u 


J 


210.0 


1 fi 4- 
-LU.UOI zc. 





030 


QPCft Q J_ Q 7 

oou . zc y . 1 


Lick 


2714.0 


K' 


210.0 


14.143 ± 





100 x 


1 ^03 zL 4- 1 39 3 X 

IOUO.1 ZC IOZ.O 


Lick 


2997.0 


K' 


210.0 


14.321 ± 





100 x 


1 97ft 1 4.IIO QX 
1Z 1 U. 1 ZIZ lli.O 


Lick 


3279.0 


K' 


210.0 


14.384 z\z 





100 x 


1 904 1 4- 1 Ofi X 

liU4 . X ZIZ XIJU.U 


Lick 


3^9 
ouuz . u 


K' 


210.0 


14.518 ± 





100 x 


1 Ofizt 3 4- Q3 fi x 
iuui. zc yo.u 


Lick 


3846.0 


K' 


210.0 


14.583 ± 





100 x 


1 009 ^ 4- 8S 9 X 


Lick 


4129.0 


K' 


210.0 


14.650 ± 


0. 


100 x 


Q49 PC 4_ 09 QX 

yiz.o zc oz.y 


Lick 


4413.0 


K' 


210.0 


14.818 ± 





100 x 


Q07 4 4- 71 O x 

OU 1 . 1 ZC I l.U 


1V1 dtillLlXIl 


1 090fi n 

XUZrUU. u 


j 


600.0 


17.760 ± 





059 


1 Q9 4 4_ 7 

IOZ . 1 ZC f . u 


A/In o~Tn i tti 

1V± (Xtill Li 111 


1 0396 n 


R 


600.0 


91 sp;n 4- 





390 


6 69S 4- 9 000 


l\/T Q rrn 1 1 vvi 


1 i ^9^ n 

X lUiU. u 


Y 


300.0 


1 S 40K 4- 





183 


Qft 4Q 4_ 1 4 Qft 

yu.10 zc n.yu 


M 11 um 


11526.0 


I 


600.0 


19.880 ± 





.140 


30.49 ± 3.69 


Magnum 


12846.0 


K 


480.0 


16.491 ± 





086 


172.9 ± 13.2 


Magnum 


12846.0 


R 


600.0 


> 21.28 






< 11.20 


Magnum 


14106.0 


H 


540.0 


17.597 ± 


0. 


082 


96.87 ± 7.05 


Magnum 


14106.0 


I 


600.0 


20.220 ± 





200 


22.29 ±3.75 


Magnum 


15366.0 


J 


540.0 


18.320 ± 





103 


79.08 ±7.16 


Magnum 


15366.0 


R 


600.0 


> 21.36 






< 10.41 


Magnum 


16506.0 


Y 


540.0 


19.750 ± 


0. 


430 


28.02 ±9.16 


Magnum 


16506.0 


I 


600.0 


20.460 ± 





280 


17.87 ±4.06 


Kuipcr 


5098.5 


I 


1055.0 


18.452 ± 





085 


113.6 ±8.6 


Kuiper 


11260.0 


I 


2176.0 


19.798 ± 





010 


32.88 ±0.30 


Kuiper 


15382.0 


I 


3530.0 


> 19.65 






< 37.68 


Kuipcr 


18584.5 


I 


2701.0 


> 19.23 






< 55.48 
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Table 5. PAIRITEL JHK S Secondary Standards 



a S J cr H cr K a cr 

deg deg mag mag 



355.107649 


31.795298 


11.681 


0.012 


11.445 


0.011 


11.418 


0.007 


355.058508 


31.790998 


12.734 


0.007 


12.426 


0.007 


12.385 


0.011 


355.037846 


31.737404 


13.095 


0.007 


12.717 


0.007 


12.627 


0.010 


355.037722 


31.708279 


13.207 


0.052 


12.952 


0.048 


12.872 


0.036 


355.134554 


31.744791 


13.704 


0.031 


13.393 


0.029 


13.343 


0.024 


355.077212 


31.725382 


13.832 


0.024 


13.567 


0.014 


13.506 


0.022 


355.105243 


31.771235 


14.137 


0.007 


13.762 


0.007 


13.709 


0.015 



Magnitudes are observed values, not corrected for Galactic 
extinction, and are reported in the 2MASS (Vega) system. 



Table 6. Optical-IR Secondary Standards 



a 




8 




R 


I 


Y 




J 


H 


K s 


deg 




deg 


mag 


mag 


mag 


mag 


mag 


mag 


355. 


016646 


31 


770287 


17.046 


16.379 


16 


.095 


15.678 


15.203 


15.069 


355. 


.025377 


31 


744127 


13.699 


12.469 


11 


.816 


11.297 


10.683 


10.476 


355. 


.034260 


31 


.771452 


18.039 


16.210 


15 


.318 


14.824 


14.211 


13.949 


355. 


.037846 


31. 


.737404 


14.231 


13.684 


13 


.418 


13.095 


12.717 


12.627 


355. 


.049875 


31 


.745661 


15.978 


15.439 


15 


.211 


14.942 


14.622 


14.563 


355. 


058508 


31 


790998 


13.740 


13.229 


12 


.990 


12.734 


12.426 


12.385 


355. 


.058815 


31 


780569 


16.792 


16.262 


16 


.043 


15.775 


15.457 


15.385 


355. 


.066002 


31 


.793428 




15.960 


15 


.668 


15.310 


14.894 


14.839 


355. 


.068623 


31. 


725666 


17.911 


16.674 


17 


.607 


16.794 


16.042 


15.371 


355. 


.077212 


31 


.725382 


14.820 


14.314 


14 


.042 


13.832 


13.567 


13.506 


355. 


105243 


31 


.771235 


15.358 


14.765 


14 


.452 


14.137 


13.762 


13.709 



Magnitudes are observed values, not corrected for Galactic 
extinction, and are reported in the Vega system. 



